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We investigate the transfer of entanglement from source particles (SP) to target particles (TP) in 
the Heisenberg interaction H = si • S2. In our research, TP are two qubits and SP are two qubits 
or qutrits. When TP are two qubits, we find that no matter what state the TP is initially prepared 
in, at the specific time t — tt, the entanglement of TP can attain to 1 after interaction with SP 
which stay on the maximally entangled state. For the TP are two qutrits, we find that the maximal 
entanglement of TP after interaction is relative to the initial state of TP and always cannot attain 
to 1 to almost all of initial states of TP. But we discuss an iterated operation which can make the 
TP to the maximal entangled state. 
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INTRODUCTION 



Entanglement has been considered as a physical re- 
source of quantum information processing. It is pro- 
foundly important in quantum teleportation [1], quan- 
tum computing [2], cryptography [3], and quantum 
games [4,5]. Using Jayness-Cumming-type interaction 
between atoms and cavity fields [6] , entanglement trans- 
fer from two-mode squeezed vacuum state to two sepa- 
rable atoms which are prepared in pure state has been 
studied in Ref. [7]. Recently, Zou et al have discussed 
the entanglement transfer from some entangled two-mode 
fields to mixed qubits [8] and Zhang et al have discussed 
the entanglement transfer from photons to atoms [9] . 

In the paper [9], they conclude that the maximally 
entangled state of photons can lead the atoms to the 
maximal state of entanglement by nonlinear interaction. 
That is, the entanglement among photons transfers to the 
atoms by the interaction. Similarly, the transfer process 
between atoms and atoms stimulated our interests and 
then we research it. Fortunately, we find some interesting 
things in the cases of two pairs of qutrits. In the process 
of study, we consider the Heisenberg interaction — H = 
§i • S2, and mainly focus on the entanglement transfer 
from source particles (SP) which are two qubits or qutrits 
to the target particles (TP) which contain two qubits. 
When SP are two qubits, the process of the transfer is 
clear and easy to understand, because the quantity of 
entanglement between two qubits is well-defined. But 
when SP are two qutrits, we have not a good definition 
of entanglement degree and only have two invariants I\ 
and I2 [10] . In the second section, we introduce the model 
and the study method. In the third section, the transfer 
between two pairs of qubits is studied. And in the Fourth 
section, we study the cases of transfer from two qutrits 
to two qubits. The Final section is the conclusion. 
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FIG. 1: The picture of the Model. 



THE MODEL AND THE RESEARCH METHOD 

As shown in Fig.l, there are four particles in this 
model, and two of them 1 and 2 are called target par- 
ticles (TP) which are two qubits and initially stay on the 
pure state \1p12} 
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111) 



(1) 



The other two particles 3 and 4 are called source parti- 
cles (SP) which are two qubits or two qutrits and initially 
stay on the pure state \1p34). The initial density operator 
of the whole system is 



p(0) -Pi 2 (0)®p 34 (0) = hMftM ® IV^WmI (2) 

where pi2(0) is the initial density operator of TP and 
^34(0) is similar. Then, particles 1 and 3 interact under 
the interaction Hamiltonian H — si ■ S2 and particles 2 
and 4 are also acted with H. 
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The time evolution operator of the whole system can 
be derived U(t) — u\ 3 {t) (g> u 2 i{t), where u\ 3 {t) is the 
time operator between particles 1 and 3, then u 2 i(t) is 
similar. Now we can obtain the density operator p(t) of 
the system at time t 

p(t) = U(t) P {Q)U\t) 

= u 13 (t) <g> u 2i (t)p 12 (0) <g> p 3i {0)u\ 3 (t) <g> ul 4 {t\3) 

To get the entanglement between two target particles, 
one needs to trace over the SP 3 and 4 to obtain the 
reduced density operator for the TP p\ 2 {t) 



Pi2(t) = Tr 3A p(t) 



(4) 



From (1), we can calculate the initial quantity of en- 
tanglement £i 2 (0) of TP 



E 12 (0) = 2 1 sin 6>i cos0i| 



(5) 



where we use the Negativity [11,12] to denote the quan- 
tity of entanglement. For a bipartite system described 
by the density matrix, the negativity criterion for entan- 
glement of the two subsystems is given by the following 
formula: e — — 2 J2i ^7 where the sum is taken over the 
negative eigenvalues of the partial transposition of the 
density matrix p. The value e = indicates that the two 
subsystems are separable. This function varies between 
and 1, and monotonically increases as the entanglement 
grows. Then from (4), we can calculate the Negativity 
of TP at the time t — E 12 {t). We can know the detailed 
process of the entanglement transfer from SP to TP by 
studying the E\ 2 (t). 



ENTANGLEMENT TRANSFER FROM TWO 
QUBITS TO TWO QUBITS 

The initial state of SP IV'34) is 

1^34) = cos0 2 |OO)+sin0 2 |ll) (6) 
And the interaction Hamiltonian between 1 and 3 is 



H = si ■ s 3 = 



( 1 \ 
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The time evolution operator between particle 1 and 3 
u 13 (t) is 



Ul 3 (t) 



(I 

\{l + e lt ) 

-i(-l + e rf ) 
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The interaction Hamiltonian and the time evolution 
operator between particle 2 and 4 are similar with parti- 
cle 1 and 3. 

The initial quantity of entanglement of SP £34(0) is 



E 34 (0) = 2 1 sin 2 cos (9 2 1 



(9) 



After some calculations, the reduced density operator 
P\ 2 {t) for TP is derived 



Pl2{t) 



( A{t) F(t) \ 

B(t) 

C(t) 

\ F*(t) D(t) ) 



(10) 



where the matrix basis is chosen as 
{|00),|01),|10),|11)}. The coefficients in Eq. (10) 
are functions of time and given by 



A{t) = (cos 2 sin 0i sin 2 ^ — sin 2 cos 0i cos 2 ^) 2 

+ COS 2 01 COS 2 2 

B(t) = - sin 2 t sin 2 (0i + 2 ) 
C{t) = B{t) 

D{t) = (cos 02 sin 0i cos 2 -— sin 2 cos 0i sin 2 -) 2 

+ sin 2 0i sin 2 2 
F(t) = e- Jt [-cos02sin0 2 sin 2 ^(cos 2 0i +e 24t sin 2 0i) 

+ cos 0i sin 0i cos 2 ^(cos 2 2 + e 2U sin 2 2 )] (11) 

And the analytic expression of the quantity of en- 
tanglement of TP at the time t, E\ 2 (t) is very com- 
plicated. And one who wants to know its real form, 
can calculate it using the following formula — Ei 2 (t) — 
Max {0, y/(B -C) 2 +4|F| 2 - (B + C)}. From now on, 
we use the numerical method to study the problem and 
give some numerical results. 

First, we find that Ei 2 (t) has a period that is 2-7T, that 
is to say 

E 12 (2kir) = E 12 (0) = 2|sin0icos0i| (12) 

where k is the natural number. 
Second, when t = it, 

£12(71-) =£34(0) = 2| sin 2 cos 2 1 (13) 

That is to say, the entanglement of TP can attain to the 
entanglement of SP. Assuming that the entanglement of 
SP is 1 and initially the entanglement of TP is very small, 
we can enhance the entanglement of TP and make them 
stay on the maximal entangled state under the interac- 
tion H. 
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FIG. 2: The entanglement of 1,2 particles S12 as a function of 



t. (a).0i =0,02 = f ; (b).0i 
(d).0i = f,0 2 = 



(c). 



Finally, we give some numerical pictures to show the 
propriety of the cases. 

From the Fig.2, one can see that when £12(0) = 
0,-^34(0) = 1 and as time goes by, the entanglement of 
TP increases from to 1 when t = 7r; then decreases 
from 1 to 0. To the picture (b), (c) and (d), the initial 
entanglement £12(0) of TP is not equal to 0. The process 
of transfer contains five sections, (i). the entanglement 
of TP E12 decreases from £12(0) to zero; (ii). E± 2 stays 
on the zero during a period of time; (hi). £12 increases 
from zero to £34(0) when t = ir and then decreases from 
£34(0) to zero; (iv). £12 also stays on the zero during 
a period of time; (v). Ei 2 increases from zero to £12(0) 
when t = 2n. 



ENTANGLEMENT TRANSFER FROM TWO 
QUTRITS TO TWO QUBITS 

In this section, the SP are two qutrits and so we re- 
search the transfer process from two qutrits to two qubits. 
The initial state of SP IV'34) is 



1^34) =fc |00)+fci|ll)+fc 2 |22) 



(14) 



where kg,ki, k 2 satisfy the normalizing condition |fco| 2 - 

N 2 + M 2 = i. 

And the interaction Hamiltonian between 1 and 3 is 



H = Sf s 3 = 



(* 
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The time evolution operator between particle 1 and 3 
u 13 (t) is 



Ul 3 {t) = 



/ii 0\ 

x 2 x 3 

Xi x 3 

x 3 X4. 
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x 2 
Xl j 



(16) 



where x\ 



3 , %i — 3 

Just as the last section, particle 2 and 4 do the same 
operation with particle 1 and 3 in the interaction Hamil- 
tonian and time evolution operator. The two invariants 
of IV34) is [10] 
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After some calculations, the reduced density operator 
Pi 2 {t) for TP is derived 



/ a(t) e(t) \ 

b(t) 

c(t) 

V e*(t) d{t) J 



(18) 



where the matrix basis is chosen as 
{|00),|01),|10),|11)}. The coefficients in Eq. (18) 
are functions of time and given by 
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a(t) 

b(t) 

c(t) 
d(t) 

e(t) 



Cg COS 2 ' 

~— 3it 



-3 it 



fc^j cos 2 6»i + e —— [(2 + ) 2 ki cos 6>i + 2(-l + j 2 k sin 6>i] [(1 + 2e^ ) 2 fc x cos 6>i + 2(-l + e" 2 " ) 2 /c sin 6>i] + 
81 



e -■■ r ,„ 3n. 
~81 
2e" 3it 



81 

2e- 3it 



[(2 + e^) 2 fc 2 cos 6»i + 2(-l + eT") 2 ^ sin6>i][(l + 2e 1 ^) 2 fc 2 cos 6^ + 2(-l + e^) 2 k 1 sin6q] 

(-1 + e 2 5 i ) 2 [(l + 2e 3 s i ) 2 fci cos6»i + (2 + e^) 2 fc sin^][(2 + e^fk x cos6>i + (1 + 2e^) 2 k sin6>i] 
(-1 + eT) 2 [(l + 2eT i ) 2 fc 2 cos6>i + (2 + e^) 2 /^ sin6>i][(2 + e^) 2 k 2 cosd 1 + (1 + 2e^) 2 k sin0i] 



81 

6(i) 

fc| sin 2 6»i + ^— — [2(-l + e 3 5 i ) 2 fci cos6»i + (2 + e^ L ) 2 k sin0i][2(-l + e^) 2 k x cos6>i + (1 + 2e^r) 2 k sin0i] + 
'81 



-3it 



81 



- [2(-l + ) 2 fc 2 cos 6»i + (2 + e^fh sin 0i] [2(-l + ) 2 fc 2 cos X + (1 + 2^ ) 2 fci sin 0i] 
-fc cos0 x [2(-l + e 3 s i ) 2 fci cos 6»i + (2 + e^) 2 fc sm0i] 



fc 2 sin0i[(l + 2e 2 ) fc 2 cos0i +2(-l + e 2 ) fcisin0i 



-3it 



H — - — [(2 + e~) 2 fci cos0i + 2(-l + e^) 2 k sin0i][2(-l + e~) 2 k 2 cos0i + (1 + 2e~) 2 fci sin0i] 
81 



(19) 



Compared with the last section, the analytic expression 
of the quantity of entanglement of TP at the time t — 
E\ 2 (t) is more complicated. We still use the numerical 
method to study the problem and give some numerical 
results from now on. 



First, we find that E 12 (t) has a period that is ^ 
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= E 12 (0) = 2 1 sin 0i cos 6<i | 



(20) 



where k is the natural number. 

Second, being similar to the case that SP are two 
qubits, when t = T/2 i.e. t — 2ir/3, SP can influence the 
entanglement of TP to a great extent. So we can calculate 
the analytic expression of the reduced density operator 
Pi 2 (t) and the entanglement at the time t = ^p, denoted 
as pi 2 ( 2 | : ), Ei 2 (^-). Obviously, Si 2 (^ r ) is a function of 
{0i, fco, k\, k 2 } and by some transformations the parame- 
ters can become {0i , I\ , I 2 }. For a given value of 0i , there 
is a pair of I\,I 2 that makes the entanglement Ei 2 ( 2 ^-) 
to the maximal quantity. Considering the complexity of 
analytical form of £a 2 (^), we have to use the numerical 
method to find the maximal quantity of Ei 2 (^-), denoted 
as E 12Max (?f), and 7 1 ,J 2 matched with E 12Max (^f) to 
a special 6\. 

Now, we introduce how to find the Ei2Max\^f)- 
Firstly, the area to satisfy the physical condition is very 
small in the picture of number pair I\ — I 2 . To show it 
clearly, we need to make a transformation 



i, 



h 

h - -h 
2 



(21) 




FIG. 3: 
with 6i 



The picture of I 2 — Ii and 9 maximal point matched 
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The physical region transform in the Fig. 3 and A, B, C 
three points are: 
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From Fig. 3 one can see the 9 maximal points are 
almost on the line AB. The real equation of frontier AB 
is very complicated, so we give an approximate one. That 
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is 

4 = - 1 (22) 

2 3 1 6 v ; 

And the error is small or equal to 1%. So we can 
infer that to a special 6\, the point — {/i,^} that makes 
-^i2(^f ) attain to the maximal quantity is approximately 
on the line AB. Now, you may ask that how Q\ is matched 
with I\,l2- By some calculations, we can get a better 
approximate result 

Ii = -0.08756 sin 2 26>i - 0.07911 sin 26»i +0.5 (23) 

Its error is small or equal to 1.4%. Now, if you give a #i, 
you can first use Eq.(23) to get I\ and then lead the re- 
sult you calculated to the Eq.(22) and (21), then you get 
the ii,/2 which matched with the maximal EyiMax^f- )■ 
And the value of E\iMax(^) can be obtained. 

Secondly, from the above table, we know point A and 
B are very special in all of the states. State A is the 
maximal entangled state of two qutrits. And the form of 
state B is likely with the maximal entangled state of two 
qubits. But the state B cannot elevate the entanglement 
of TP to unity and only to 0.8. 

Now, some interesting things happened. To the case 
that SP are two qubits, the entanglement of TP after in- 
teraction must attain to 1 when SP are in the maximally 
entangled state. But the case that SP are two qutrits is 
very different from two qubits. The maximally entangled 
state of the qutrits cannot make all states of TP attain 
to the maximal state as the Fig. 4. shown. 

More interesting is that there does not exist a state 
that can make the entanglement of state of TP to unity 
in all states of qutrits. That means that there is a special 
state matched with a special 6\ that can make the entan- 
glement of states of TP to a maximal quantity which less 
than 1. The details have been shown in Fig. 4. 

From the Fig. 4, we can see that all states of qutrits 
cannot elevate the entanglement of TP to 1 except the 
case that when SP is on the state A and the initial state 
of TP is maximally entangled state. But, if you want 
to make the entanglement of TP attain to 1 only us- 
ing state A, you can make many times just like the fold 
line in Fig. 4 has shown. It is after interacting with SP 
which stay on state A that the entanglement of TP is 
from initial 0.2 to 0.62. Similarly, with 0.62 as the initial 
value, the entanglement of TP is from 0.62 to 0.82 after 
interacting with SP which stay on state A again. Doing 
this operation iteratively, after interacting with SP four 
times, the entanglement of TP can attain to 0.96. 

DISCUSSION AND CONCLUSIONS 

In the interaction Hamiltonian H = s% • s%, the entan- 
glement transfer from SP to TP is investiged. We mainly 
focus on two cases. 
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FIG. 4: Black real line denotes the initial entanglement of 
TP; dashed line denotes the entanglement of TP when i = ^ 
after interaction with SP which stay on the A state initially; 
dotted line denotes the entanglement of TP when t = 4p after 
interaction with SP which stay on the B state initially; and 
the red real line denotes the maximal entanglement quantity 
of TP in all states of qutrits after interaction; the fold line 
denotes the process that entanglement of TP initially is 0.2 
and attains to 0.96 after interaction with SP 4 times. 

Case 1: When SP are two qubits, we can enhance the 
entanglement of TP to 1 by interacting with SP whose 
entanglement is 1. 

Case 2: When SP are two qutrits, generally, we can- 
not enhance the entanglement of TP to 1 by the interac- 
tion with SP, even if SP is on the maximally entangled 
state. More interesting thing is that to a special initial 
entanglement of TP, there exists a state in all states of 
qutrits which satisfy the physical condition. Compared 
with other states of qutrits, this state can make a largest 
enhancement to the entanglement of TP after interaction. 
Generally, this state is not the maximally entangled state, 
but it can make TP to a larger entangled state than the 
maximally entangled state did. 

Considering this thing carefully, we find that in the 
case 1, the maximal entanglement of TP after interaction 
has no relation with the initial state of TP and always 
can attain to 1 to all initial states of TP, as long as SP 
stay on the maximally entangled state initially. In the 
case 2, the maximal entanglement of TP after interaction 
is relative to the initial state of TP and always cannot 
attain to 1 to all states of TP except one condition that 
the initial state of TP and SP is the maximally entangled 
states. But this case is useless. Because the initial state 
of TP is already the maximally entangled state, it is not 
necessary to interact with SP. Therefore, we only have 
interest on the initial state of TP which is not maximally 
entangled state. 

From the above discussion, we find that there are some 
differences between case 1 and case 2. Firstly, the re- 
lation between the maximal entanglement of TP after 
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interaction and the entanglement of TP initially is differ- 
ent. Secondly, the behavior and function of maximally 
entangled state is different in the entanglement transfer 
process. Our next work is to give the reason that brings 
about those differences exactly. 

Now, we consider the reason from two aspects. On the 
one hand, the difference derives from the form of interac- 
tion. So we will use more general form of interaction, for 
example H — Jisfs^ + + ^3 s i s 2 + ^i s i + ^2 S I- 

And we will adjust the quantity of those parameters and 
see whether the difference will vanish or not. On the 
other hand, we conjecture that the difference of entan- 
glement between qubits and qutrits lead to those differ- 
ences. As we all know, there does not exist a well concept 
to describe the entanglement of two qutrits. We expect 
to hnd a proper and general definition of entanglement of 
two qutrits based on the entanglement transfer to obtain 
the essential reason. 

This work is supported in part by NSF of China Grant 
Nos. f0575053 and 106050f3 , Program for New Cen- 
tury Excellent Talents in University, and the Project- 
sponsored by SRF for ROCS, SEM. 
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